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Abstract 

We develop a systematic theory of multi-particle excitations in strongly interacting Fermi sys- 
tems. Our work is the generalization of the time-honored work by Jackson, Feenberg, and Campbell 
for bosons, that provides, in its most advanced implementation, quantitative predictions for the 
dynamic structure function in the whole experimentally accessible energy /momentum regime. Our 
view is that the same physical effects - namely fluctuations of the wave function at an atomic 
length scale - are responsible for the correct energetics of the excitations in both Bose and Fermi 
fluids. Besides a comprehensive derivation of the fermion version of the theory and discussion of the 
approximations made, we present results for homogeneous ^He and electrons in three dimensions. 
We find indeed a significant lowering of the zero sound mode in ^He and a broadening of the col- 
lective mode due to the coupling to particle-hole excitations in good agreement with experiments. 
The most visible effect in electronic systems is the appearance of a "double-plasmon" excitation. 

PACS numbers: 67.30.-n, 67.30.em, 71.10.Ca, 71.15.Qe, 71.45. Gm 



1 



I. INTRODUCTION 



This paper is concerned with a systematic theory of multi-particle excitations in Fermi 
systems. We utilize an equations of motion method that has been used in the past as a 
vehicle for many purposes: the derivation of the time-dependent Hartree-Fock (TDHF) 
theory Jlj-Q, its analog for strongly interacting systems |], and for studying single- and 
multi-particle correlations in strongly interacting Bose liquids ^. 

The simplest way to deal with excitations is to assume that the low-lying excited states 
of a quantum fluid can be characterized by the quantum numbers of a single particle. This 
is the core idea of Landau's quasiparticle picture of "normal" quantum fluids ^ as well 
as of Feynman's theory of collective modes in the helium liquids . It is appropriate for 



many long wavelengths excitations such as sound waves in Bose fluids or plasmons in an 
electron liquid. 

Already Feynman realized that this concept is insufficient to describe higher-lying excita- 
tions, most prominently the "roton" in ^He. Intuitively appealing, he introduced "backflow" 
correlations ITT]. These are recognizable as a new type of excitations, depending on two par- 
ticles: pair fluctuations . The notion is plausible: For excitations at wavelengths comparable 
to the interparticle distance, the time-dependence of a system's short-ranged structure is 
expected to be relevant. 

The presently state-of-the-art theory for Bose liquids originates from pioneering studies 
by Jackson, Feenberg p!2|-p!6[] , and Campbell and collaborators . Recently, a complete 



solution of the pair equation of motion has been accomplished in *^He [0], showing that 
the "uniform limit approximation" of Refs. |], p!2|-p^ is surprisingly good. Consequently, 
theoretical improvement must be sought in three-body and higher-order fluctuations [|18| . 

Although quite successful for bosons, there exists to-date no fermion version of the theory. 
We therefore develop here the generalization of the equation of motion method for pair 
fluctuations to fermions. We calculate the fermionic density-density response function 
r'; t — t'), relating the induced density fluctuation 6p{r;t) to a weak external perturbation 
^ext(r; t). In a homogeneous system this is written in momentum space as 

5p(q; w) = px(?; ^ext(q; , (i.i) 
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where p is the particle number per volume VL. We choose Fourier transforms 

/(r-^) = (1-2) 

q 

to have the same dimension in q- and in r-space. 

The imaginary part of xil'i ^) is the experimentally accessible dynamic structure factor, 

h 

S{q] w) = '^m[x{q; u)] 9{u) . (1.3) 

TT 

The dynamic structure factor satisfies, amongst others, the sum rules 

POO 

■mo = S{q) = / dhu S{q;u), (1.4) 
Jo 

mi = — — = / dhjj hu S {q; u) , (1.5) 
2m Jq 

where S{q) is the static structure factor. 

We develop our theory with the following objectives: 

• Technically, the extension of the Jackson-Feenberg-Campbell theory to Fermi systems 
amounts to including time-dependent two-particle-two-hole excitations. We require 
that the fermionic x{q'i ^) reduces to that of the boson theory in the appropriate limit. 

• For bosons, neglecting pair- and higher order fluctuations yields the famous Bijl- 



Feynman spectrum [10 



Its fermionic counterpart is the random-phase approximation (RPA), formulated in 
terms of effective interactions [jl9|. We require that our theory reduces to the RPA if 
pair fluctuations are ignored. This implies, in particular, that we obtain in this case a 
response function of the form 

x(g; oj) = . (1.7) 

1 - Vp-hiq) xo{q; ^) 

Here, Xo{q'i^) is the Lindhard function and Vp-h(9) an appropriately defined static 
"particle-hole interaction" or "pseudo-potential". 

One of the tasks of microscopic many-body theory is to justify and calculate effective 
interactions such as Vp-J^q), as far as this is possible. Using Jastrow-Feenberg correlation 



functions [|T3[ to tame the microscopic hard-core repulsion, it has been shown under 
what assumptions a density response function of the RPA form ( p..7| ) can be obtained, and 
a microscopic expression for the static effective interaction Vp_h(g) was derived. Under what 
conditions a form ( |1.7D is meaningful at all will be discussed in depth below. 

A phenomenological approach to define a particle-hole interaction or "pseudo-potential" 
for ^He and electrons was introduced by Aldrich, Iwamoto, and Pines [^, They deter- 
mined the physically intuitive and necessary requirements for Vp-^lq), postulating that the 



dynamic response is given by the RPA form (|1.7|) . Reflecting the same physics, the Vp-J^q) 
derived from microscopic many-body theory is very similar to the Aldrich-Iwamoto-Pines 
pseudopotentials. The microscopic derivation leads to a Vp-h(q') that is uniquely determined 
from the static structure function by the two sum rules ( |1.4| )-( pT5D . Defining the RPA this 
way leads for bosons to the Feynman approximation ( |1.6|) for the spectrum of collective 
excitations. From here on, we will use the term "RPA" and "Feynman spectrum" in this 
sense. 

Our work is organized as follows: Section || introduces the basic quantities and the most 
important tools of variational and correlated basis function (CBF) theory. For details, the 
reader is referred to review articles ||21]] and pedagogical material a brief outline of 
our notations and definitions is given in appendix ^ Section |T| is the core of our work; 
it provides the derivation of the equations of motion, including pair fluctuations. We show 
that the theory can be mapped onto a set of TDHF equations with energy- dependent, 
effective interactions. Thus, our work provides the logical generalization of Ref. |^, where 
single-particle fluctuations led to a TDHF theory with static effective interactions. 

Section |^ focuses on the practical implementation of our theory. We formulate, among 
others, the "convolution approximation" for fermions. In Section we derive the density- 
density response function x(g; u) and discuss its features. 

Modern techniques of many-body theory are robust against the details of the interparti- 
cle interaction. We can therefore use the methods developed here to examine the dynamics 
of two very different systems: The very strongly interacting ^He whose interaction is char- 
acterized by a repulsive hard core and a short-ranged attraction, and electrons with their 
rather tame but long-ranged Coulomb interaction. Section [V^ implements our method for 
bulk ^He and the electron liquid. In '^He, we compare with neutron scattering experiments 
carried out at the Institut Laue Langevin (ILL) in the group led by R. Scherm p5|^5[. The 
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energetics of the collective mode as well as the width of the spectrum at high momentum 
transfers are significantly improved compared to RPA predictions. In the homogeneous elec- 
tron liquid the pair-excitation theory predicts plasmon damping as well as double-plasmon 
excitations. Experimental verification of the double-plasmon excitation in recent inelastic 
X-ray scattering measurements has added new interest in studying the dynamics of 

electrons. 



Our results are summarized in Sec. [VI1| where we also discuss the directions of future 
work. 

Appendices give further details on the derivations, and Appendix ^ a very brief 
summary of the minimal implementation of our theory. 

II. THEORY FOR STRONGLY INTERACTING FERMIONS 
A. Variational theory 

Microscopic many-body theory starts with a phenomenological Hamiltonian for N inter- 
acting fermions, 

^ = -E|i^^' + E-(|r^-r,i) . (2.1) 

i i<j 



For strong interactions, CBF theory |T3| has proved to be an efficient and accurate method 
for obtaining ground state properties. It starts with a variational wave function of the form 

where $o(l; A^) is a model state, normally a Slater-determinant, and "i" is short 
for both spatial and u discrete (spin and/or isospin) degrees of freedom. The correlation 
operator F{1, . . . , N) is suitably chosen to describe the important features of the interacting 



system. Most practical and highly successful is the Jastrow-Feenberg |0 form 

(2.3) 



Fil,...,N) = exp{l- 
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^ M2(ri,rj)+ ^ M3(ri,rj,rfc) + . . . 

-l<j<i<A^ l<i<j<k<N 

The Un(ri, . . . , r„) are made unique by requiring them to vanish for |rj — rj| — )-oo ("cluster 
property"). 



From the wave function (|2.2[), (|2.3|), the energy expectation value 



H^,o={^o\H\^^) (2.4) 



can be calculated either by simulation or by integral equation methods. The hierarchy 
of Fermi-Hypernetted-Chain (FHNC) approximations is compatible with the optimization 
problem, i.e. with determining the optimal correlation functions M„(ri, . . . , r„) through 
functionally minimizing the energy 

-0. (2.5) 



Sun{ri, . . . ,r„) 

Due to the multitude of exchange diagrams, the Fermi-HNC (FHNC) and corresponding 
Euler equations can be quite complicated [^]; the simplest approximation of the Euler 
equations (|2.5|) that contains the important physics is spelled out in App. |A 1 . 



The optimization of the correlations also facilitates making connections with other types 
of many-body theories, such as Feynman-diagram based expansions and summations 

B. Correlated Basis Functions 

Although quite successful in predicting ground state properties of strongly interacting 
systems, the Jastrow-Feenberg form (|2.3| ) of the correlation operator F has some deficiencies. 
The most obvious problem is that the nodes of the wave function (|2.2|) are identical to those 
of the model state |$o)- To improve upon the description of physics, CBF theory 1^1], |28 



uses the correlation operator F to generate a complete set of correlated and normalized A^- 
particle basis states through 



where the {|$m)} form a complete basis of model states. Although the \^m) are not orthog- 



onal, perturbation theory can be formulated in terms of these states |13|, ^ . We review here 



this method only very briefly, details may be found in Refs. ^ and ^ the diagrammatic 
construction of the relevant ingredients is given in Ref. ^ 

For economy of notation, we introduce a "second-quantized" formulation of the corre- 
lated states. The Jastrow-Feenberg correlation operator in (|2.3| ) explicitly depends on the 
particle number, i.e. F = -FJv(l, . . . , A^) (whenever unambiguous, we omit the corresponding 
subscript). Starting from the conventional a],, a^., creation and annihilation operators a;|., 
of correlated states are defined by their action on the basis states: 

\ai ^ F,^A 1$^) / {^J\a,F;^F^^^ai\^^f' , (2.7) 
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W ^m) = F,_,a^ 1$^) / {^m\alFl^F^_^a,\^^f^ . (2.8) 

According to these definitions, a\ and a^. obey the same (anti-) commutation rules as their 
uncorrelated cousins, but they are not Hermitian conjugates. If \^m) is an A^-particle state, 
then the state in Eq. ( p.7|) must carry an (A^+l)-particle correlation operator, while that in 
Eq. ( |2.8| ) must be formed with an (A^ — l)-particle correlation operator. 

In general, we label "hole" states, which are occupied in |$o), by h, h', hi and 
unoccupied "particle" states by p, p', pi , etc. To display the particle-hole pairs explicitly, 
we will use alternatively to l^&m) the notation |\I'pi...p^/ii...h,^). A basis state with d particle- 
hole pairs is then 

= aJr--aLaft,---ahJ^o)- (2.9) 

The execution of the theory needs the matrix elements of the Hamiltonian, the unit 
operator, and the density operator. Key quantities are diagonal and off-diagonal matrix 
elements of unity and H'= H — Hq o 

= (^„,|^„) = 5„,,n + AT^.n , (2.10) 

^m,n = W^m,n + \ {H^,m + i^„,n - 2ifo,o) A^m.n ■ (2.11) 



Eq. (|2.11 ) defines a natural decomposition [31, 32| of the matrix elements of H'^^. 



The ratios of normalization integrals, /m,m = {^mW^ F\(^rn) , define the factors 



^p\...p^h\...h^ — — \J Im,m/ Io,o ■ (2.12) 

For large particle numbers and d'^N these factorize as 

= ^2llll3!± + o(iv-i) . (2.13) 

Likewise, to leading order in the particle number, the diagonal matrix elements of H' = 
H — Ho^o become additive, so that for the above d-pair state we can define the CBF single 
particle energies 

d 

(^mli^l^m) = 5^ep,,, + 0{N-') , (2.14) 

i=l 

with Cph = ep- Ch. 

For the off-diagonal elements Om.n of an operator O (specifically the Hamiltonian, the 
unit-, density- and current-operator) we sort the quantum numbers mi and Ut such that 



|\I/m) is mapped onto |\l/n) by 

l^m) = ttJnia|„2 ■■■al^^a^^--- a^^a^^ |^n) • (2.15) 

From this we recognize that, to leading order in A^, any 0„i,n depends only on the difference 
between the states \^m) and |\l/n) and not on the states as a whole. Consequently, Om,n 
can be written as matrix element of a rf-body operator 

Om.n = (mi m2 . . . nid |C(1, 2, . . . rf) |ni n2 ... nd)a ■ (2.16) 



(The index a indicates antisymmetrization.) According to (|2.16|) , W^n and A^mn define 
(i— particle operators JV and W, e.(?. 

iVm.o = Np^p2...pd hj_h2...ha,0 = {PlP2 • • • | A/'(l, 2, . . . , C?) | /li/lg • • • hf)a , 

W^m,o = W^pip2...Pd/*i/^2.../id,0 = (PlP2---Pd| >V(l,2,...,rf) |/;,i/l2.../id)a • (2.17) 

Diagrammatic representations of A/'(l, 2, . . . ,d) and VV(1, 2, . . . , c?) have the same topology 



3T| . In homogeneous systems, the continuous parts of the Pi,hi are wave numbers Pi, hj; 
we abbreviate their difference as q^. The highest occupied momentum is hkp. 

An important consideration is, for our purposes, the connection between CBF matrix 
elements, the static structure function, and the optimization conditions for the ground state. 
The static structure function S{q) = ;^(\E'o|PqP-q|^o) is routinely obtained in ground state 
calculations; for some systems it is also available from experiments. We can also write S{q) 
as the weighted average of the matrix elements ( p.l7| ), 

S{q) = Sp{q) + ^^Zpp>hh'Npp'hh',o- (2.18) 

hh' 

where 5'p(g) is the static structure function of non-interacting fermions. 

Similarly, the optimization conditions ( |2.5D for the pair correlation function can, in mo- 
mentum space, be written in terms of off-diagonal matrix elements of the Hamiltonian: 

*^ 





F^H'F [pqpq, 


- Pq+q'] 


<fo) 







5M2(q,q') 

/a, |FtFl(5 \ = ^vv'hh'iipp^hh'^^ (2.19) 





F^H'F 









hh' ^ °' ' °' hh' 

i.e. the weighted average of the off-diagonal matrix elements -f^o.pp'/i/i' vanishes for optimized 
pair correlations. Both features will provide rules for systematic and consistent approxima- 
tion schemes for the operators A/'(l, 2, . . . , ci) and W(l, 2, . . . , ci). 
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III. EQUATIONS OF MOTION 



A. Excitation operator and action principle 



To formulate a theory of excited states for strongly interacting fermions we generalize 
the ansatz ( p.2| ) in analogy to the pair fluctuations theory for strongly interacting bosons 
0, 1^, |12], p!^-p!7|. We restrict ourselves here to uniform systems. The system is subjected to 
a small external perturbation 



ifext(t) = /rf\/iext(r;t)p(r) 



(3.1) 

where p(r) is the density operator. The correlated wave function for the perturbed state is 
chosen to be 



^o(t) 

m 



exp[-iHo,ot/h] ^o{t) 
1 



/l/2(t) 



exp 



Wit) 



exp 



if/t(t) exp if/(t) 



(3.2) 



with the excitation operator 



ph 



pp'hh' 



= U,it) + U2{t). 



(3.3) 



The particle-hole amplitudes 5u^j2 (^) ^Mpp/^/j/ (t) are determined by the stationarity prin- 
ciple for the action 



it) , Hr (^) , Hp' w ' (^) = / dt m 



.(2) 



(2)* 



with the Lagrangian |l|, ^, |^ 



C{t) = {^{t)\H + H,,t{t)-ih^\^{t) 

^o{t)\H' + H,^,{t)-ih^^ \^o{t) 



(3.4) 



(3.5) 



A "boson" version of the theory is recovered when the particle-hole amplitudes Su^pfl (t) and 
Su^pjlihhi{t) are restricted to local functions that depend only on the momentum transfers 



qW = p 



(') _ hW. 



B. Brillouin conditions 



To derive linear equations of motion, the Lagrangian (|3.5|) must be expanded to second 
order in the excitation operator U {t). For the procedure to be meaningful, one should require 
that the first order terms vanish. This is, in principle, a necessary condition, however, in 
practice it is not always possible to satisfy it rigorously. 

The first variation of the energy with respect to Su^p^it) and Su^pfl*{t) is 

6{m\H'\m)' 



= Kph (3-6) 

5u(i)(t)=5M(2)(t)=o 



and its complex conjugate. This term vanishes in the homogeneous liquid due to momentum 
conservation. 

(2) 

The variation with respect to SUp^ij^y leads to a similar condition 



5{^>{t)\H'\^>{t)) 



— ^Qppihh' — (3-7) 

<5M(i)(t)=5u(2)(t)=o 



and its complex conjugate. This condition is not rigorously satisfied by a Jastrow-Feenberg 
ground state. Recall, however, that the optimization condition ( p.5|) for pair correlations 
can be written in terms of off-diagonal matrix elements of H' in the form ( p.l9|) . If the 



correlation operator F is chosen optimally, i.e. satisfying Eq. ( p.5|) for all ra, the weighted 
averages of -ffo.n vanish. This shows precisely what an optimized ground state does: The 
Jastrow correlation function does not have enough flexibility to guarantee the Brillouin 
condition (PT7|), because -^o.pp'/i/i' depends non-trivially on four momenta, whereas the two- 
body Jastrow-Feenberg function depends only on the momentum transfer. Optimization 
has the effect that the Brillouin conditions are satisfied in the Fermi-sea average. 

To make progress we must assume that in the Lagrangian terms that are linear in the 
pair fluctuations are sufficiently small and can be omitted. Likewise, we also shall assume 
that the ground state wave function ( p.3|) is well enough optimized such that three- and 
four-body Brillouin conditions are satisfied. In momentum space, these are 

(vI/o|/fK,---,PqJvI/o) = 0. (3.8) 
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C. Transition density 



The quantity of primary interest is the hnear density fluctuation induced by the external 
fleld Hcxtit)- We regard this density as a complex quantity; it is understood that the physical 
density fluctuation is its real part. Assuming the excitation operator ( |3.3| ), it is 

5p{r;t) = 5^(vI/o|/)(r)-p|vl>p,)5««(t) 

ph 
1 

2 



pp'hh' 



(3.9) 

ph pp'hh' 

The matrix elements of the density, po,p/i(r) and po^pp'hh'i,'^) with respect to the correlated 
states can also be written as linear combinations of the matrix elements p^p^i^) with respect 
to uncorrelated states, and one-, two-, and three-body matrix elements of the unit operator. 
For the sake of discussion, let us briefly neglect the pair amplitudes. Since the density 
operator is local, we can commute p(r) to the right or to the left of the correlation operator 
F. The form obtained by commuting /)(r) to the left is 



Po,ph(r) = Ypl,p'h'ij^) Mp>h',ph = Po,ph(r) + YPo,p'h'{r) Np^h'^ph 

p'h' p'h' 



(3.10) 



where p^p/j(r) = Zph ($o| p(r)— p lapfl/j^o) = -Zph {h\Sp{r)\p) are, apart from the normalization 
factors Zph, the matrix elements of the density operator in a non-interacting system. 
The second form is obtained by commuting p(r) to the right of F: 

Po,ph(r) = -rPo,pft(r) + Yl ^o,pp'hh'Pl'h'fli'^) ■ (3-11) 
ph p,h> 

These two seemingly different expressions are identical, the different analytic forms appear 
only because the second quantized formulation hides the fact that the density operator is 
local. We will see below that both forms are useful. 

Including pair fluctuations, the fluctuating density (|3.9|) can generally be written as 



5p(r;t) = J^Po,p,,(r) 

ph 



y^,Mph^p'h'Su^^l,{t) + ^ Mph,p'p"h'h"Su''^l^,^,^„{t) .(3.12) 

.p'h' p'p"h'h" 

A key step that simplifles the structure of the equations of motion signiflcantly is to 
introduce a new one-body function. In analogy to the boson theory [0, we deflne new 
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particle-hole amplitudes Sv^j^{t) through 



such that 



ph 

php'h' 



(3.13) 



(3.14) 



This implies 

Yl ^ph,P'h' (t) = Y Mph,p'h' Su^jil {t) + l Yl ^ph,P'p"h'h" Sufp"h'h" (^) • (3-15) 



p'h' p'h' 

r(i) 



p'p"h'h" 



Defining M^^l^y,,^,^,, via 



^ph,p'p"h'h" — '^ph,pihi-^pihi,p'p"h'h" 
Pih\ 



(3.16) 



we can formally solve for SVpi^{t): 

SvSit) = + 1 Klp'p"H'h" ^^%H'H"it) . (3.17) 

p'pi'h'h" 

For this operation, the inverse of Mph^p'h' seems to be needed. As its calculation is not 
immediately obvious, we hasten to note that m'^^I^,^,,^,^,, is, in terms of Jastrow-Feenberg 
diagrams [0, a proper subset of the diagrams contributing to Mp/j py//i//j//. We will dis- 
cuss the diagrammatic analysis of po,p/i(r) in App. P 1| . The diagrammatic construction of 
p"h'h" spirit of Eq. ( p.l6| ) is carried out in App. 



D. The Lagrangian 



We split the Lagrangian (|3.5|) as £(t) = Cext(t) + Ct(t) + Cintit), with 



dt 



H' 



A(t) = (Mt) 



^o(t); , 



£ext(^) is obtained directly from the transition density: 



(3.18) 
(3.19) 
(3.20) 
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ph 



pp'hh' 



= ^eJ2j /iext(r; t) po,ph(r) 6v^^^{t) . 

ph 

The time- derivative term Ct{t) is, to second order in the fluctuations, 



(3.21) 



(3.22) 



+ 2E'^4VW (^'o(t)|«K'«/.'«/.^o(t)) 



+ 9 X] ^"^pXh' it)Mpp'hh',p"h"SUp,}h>' ('^) + 7 X] ^'^pp'hh' it)Mpp'hh',p''p'''h''h'''Su'j}'p'''h''h''' (t) 



Introducing the new amphtudes Sv^j^{t) defined in Eq. ( |3.13D ehminates the terms that 



couple the one- and the two-body amplitudes: 



where 



^'"ph i^)^Pf>-<P'h'^'^p'h'i^)~^ A ^'^pp'hh'i't)^pp'hh',p''p'''h''h'''^'^pi^p''' h''h''' 



(3.23) 



^pp'hh' ,p"p"'h"h"' ~ ^pp'hh',p"p"'h"h"' ^pp'hh' ,pihi^pihi,p2h2^p2h2,p"p"'h"h"' ■ (3.24) 

piP2hih2 



The second term in Eq. ( p.24| ) cancels, in a diagrammatic expansion, some terms from the 
first one (c/. App. WT). From Eqs. ( ^.21 ) and ( |3.23|) , the advantage of introducing the new 
particle-hole amplitudes Sv^i^{t) becomes obvious. 

The contributions to the interaction term are classified according to the involved n— body 
fiuctuations [/„ as defined in (lO), 



Ant(t) = /:Snt^(t) + /:|^t^(t) + /:l^^(t) 



(12), 



,(22), 



(3.25) 



with 



4n?W = ^(^o 



Ul{t)Ul{t)H' + 2Ulit)H'Ui{t) + i/'f/i(t)f/i(t)] l^o) , 
Ul{t)Ul{t)H' + Ul{t)H'U2{t) + Ul{t)H'Ui{t) + H'Ui{t)U2{t)\^ |^o) , 
Ul{t)Ul{t)H' + 2Ul{t)H'U2{t) + H'U2{t)U2{t) \ |^o) • (3.26) 
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If the Brillouin conditions ( p.6|) -( pl7D as well as their generalizations to higher order 
fluctuations were satisfled exactly, all contributions to containing [//(t)f/j(t) and 

Ui{t)Uj(t) would be zero. For fermions with optimized Jastrow-Feenberg wave functions it 
is only true in the averaged sense ( |2.19| ). These terms are nevertheless expected to be small 
in C[l^\t) since neglecting these terms is equivalent to negligible four-body correlations. 
Such a simplifying assumption is not necessary in C[l^\t) and C[ll\t) although we will see 
that the terms containing Ui{t)U2{t) and Ul{t)ul{t) in c[l^^\t) are indeed negligible. We 
keep these terms for the time being since it will turn out that their omission will suggest, 
for consistency reasons, further simpliflcations. 

The next step is to express the interaction term (|3.26|) in terms of the CBF matrix 



elements introduced on section [II B| . In the following it is understood that we sum over all 
quantum numbers when no summation subscripts are spelled out. 

^(t) = lj2Ki*(t)6u%:{t)H;p,,,,, + c.c. + \j2^^^^^^ '(3-27) 
+ ^ 5Z ^^'ph*i't)Kh,p'p''h'h''^u'^j^l''h'h''it) + c-c. , (3.28) 



Ant^^) ~ ^''^jp'hh' i't)^'^i"p'"h"h'" {t)Hpp'p"p"'hh'h"h"\0 + C-C. 

+ — ^'^lp'hh'i^)'^pp'hh',p"p"'h"h"'^'^i'^p"'h"h"'i^) ■ (3.29) 



Substituting dv^i^ (t) for du^^fl (t) leads to new coefficient functions in the interaction part of 
the Lagrangian: 

Ant(t) = Aif^(t) + <^"\t) + C.!^\t) (3.30) 

with 

Aif ^(i) = \ E (^) Kv'HH',, + c.c. + \ J2 Sv^p>:*it) Kh,p'h' ('3-31) 

-^int \^) — g Z-^^'^ph y'') ""'p'p"h'h"\'') -^^p'p"h'h",0 ' 

+ IYI ^^ph*(''^)^ph,P'p"h'h"Su^p'l"h'h"it) + c-c. (3.32) 

.'(22),. _ ^V-. (2)* (t) K^PP'^'''^ +CC 

Mnt l^J — oo / ."'pp'hh'y^J ""'p"p"'h"h"'\'') ^^p"p"'h"h"' ,0 C.C. 

+ ^^pXh'i't) Kpp>hh',p"p"'h"h"' ^u'^''p"'h"h"'i't) ■ (3.33) 
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The new coefficients Km,n are 

Kph,p'p"h'h" = H'^h,p'p"h'h" ~ '^^'ph,pihi ^pihi,p'p"h'h" ' 



Plhi 



^p'p"h'h",0 — -^pp'p"hh'h",0 Z_^-^phpihi,0 ■^^^p'p"h'h",pihi ' 

Plhi 



(3.34) 
(3.35) 



Kpp'hh',p"p"'h"h" 



Hpp'hh',p"p"'h"h"' 



2_^y^^pp'hh',pihi^pihup"p"'h"h"' ^ ^pp'hh',pihi^^^pihi,p"p"'h"h"' J 



Plh\ 



^^^pp'hh',pihi-^pihi,p2h2-^^^p2h2,p"p"'h"h'" ' 

PihiP2h2 



(3.36) 



and an analogous term for K^Fjl!l),f^iii Q ■ 



E. Equations of motion 

With the sole approximation to neglect the terms proportional to U2{t)U2{t) and 
U2{t)U2{t) , the Euler equations become 

d 



J2[^^Mph,p'h'g^-H'pf^p,^, 6v^^l{t) - J]^pp'w,o<^^^i'l'(^) 

[-^p/»yp"/»'/i"<^'"i'p"fe'?i"(0 + ^^p''h'h",o^'^pi'h'h"i^) 

- ^ [^^Mpp,^^, p,y„f^„f^,„— - Kpp>hh',p"p'"h"h'" 5u^„p,nf^ny^,„{t) 

- [Kpp'hh',p"h" Svl%„{t) + K'^,hh'fl ^^rh"it) 



= . 



(3.37) 

2 / d^r Pph,o{r)he^t{r;t) , 



(3.38) 



The time dependence of the external field can be assumed to be harmonic, with an infinite- 
simal turn-on component that determines the sign of the imaginary part 



/iext(r; t) = /iext(r; cu) [e"^* + e"''^*] e''*/^ . 



(3.39) 



This imposes the time dependence 



+ 



(3.40) 



Defining 



Epp'hh' ,p"p"'h"h"'{^) = {f'^ + ^V)^pp'hh',p"p"'h"h"' ~ ^pp'hh',p"p"'h"h"' (3-41) 
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the equations of motion for the pair fluctuations are 
1 



y^^Epp:hh',p"p"'h"h"'{ '^Wp"p'''ft"ft"'(^) - Kpp>hh',p"h" i^^p/zt" (^) + -^pp'/!h'!o^^p''/i''i^) 



Epp'hh',p"p"'h"h"' ( ^) ^'^\i"j"h"h"' (^) ~ ^^pp'hh',p"h" ^'^p"j' (^) + ^^p'hhio i'^) 



(3.42) 



All pair quantities are symmetric under the interchange of the involved pair variables, 
e.g. {pp',hh') o {p'p,h'h). We can utilize this feature to replace the fully symmetric 
Epp'hh',p"p"'h"h"'{^) by an asymmetric form, e.g. ( PT| ) which removes the factor 1/2 in Eq. 

(2±) 



The pair equations ( |3.42| ) are now solved for the ^ul ,li^,{uj) and the solutions are inserted 



into the one-body equation. The latter retains the structure of a TDHF equation, but 
with the matrix elements of H' supplemented by frequency-dependent terms. We adapt the 
definition of W^^n in ( p.ll|) by adding these corrections: 

Wph,p'h'{(^) = Wph,p'h' + ^Kph,p^p2hih2 Ep^l^^^h2,p[p'2h[h'S^)^P[p'2h^^^^ 

+ '^^p[v2hxh2fi^lj2hxh2,p[p'2Kh'Jy~^)^Xp'^^^^^ (3.43) 



+ X]-^piP2/iife2,o-^pipthi/i2,p'iPi,/i'ife^( ^)^p[p'2Kh'^,p'h' ■ (3.44) 

This TDHF form results also if the terms containing U2{t)U2{t) are retained, but the ex- 
pressions for the dynamic parts of the 11^-matrices become lengthier. 
The equations of motion for the particle-hole amplitudes are then 



2 / d^r /iext(r; uj) po,ph{r) = ^ |^(/za;+ir7) Mph^piy - 6p^p'6h,h' eph 

'' p'h' 



-E 

p'h' 
p'h' 



Wph,p'h'{^) + - (eph + Gp'h') Nph^p'h 



Wppihh'fli^) + 2 i^ph + ^p'h') Npp'hh'fi 



5v%}{uj). (3.45) 



F. Supermatrix representation 



We can now carry out exactly the same manipulations as in previous work |^ and re- 
duce these equations (|3.45 ) to the form of TDHF equations with energy-dependent effective 
interactions. 
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Equations ( p.lO|) and (|3.11|) express the density in terms of CBF matrix elements in two 
different forms. For the present purpose, it is convenient to use these two representations 
symmetrically, 



1 + 



Po,ph(r) + IY1 [Po,P'h'Np'h',ph + Po*p'h'{'^)No,pp'hh'] ■ (3.46) 



p'h' 



Using Eqs. ( |3.13|) and ( |3.40| ), the density fluctuations can then be written as 



ph 



- 2 ^ [^^'P^ 

ph 



(r) + Po^phir) Scl'f^ '{u) 



(3.47) 



(c/. ( p.lUD for the definition of Pop/i(r)). This defines new amplitudes 5c'pi^\uj). These 
relate, apart from the normalization factors, the observed density to the matrix elements of 
the density operator in the non-interacting system. The equations of motion can now be 
simplified by introducing a "supermatrix" notation. Particle-hole matrix elements together 
with their complex conjugate are combined into vectors, e.g. 



Po,ph 
P0,ph 



6c = 



6q 
5c 



(1+) 

ph 

(1-) 
ph 



(3.48) 



(and analogously for 5v^p^'^). Equation ( |3.47| ) then simply reads 



5p(r;w) = \ 5c{uj) ■ p^(r) 



The matrices 



and 



N 



Nph,p'h' Nppihh',0 
No,pp'hh' Np/h',ph 



(3.49) 



(3.50) 



^ph 







1 



1 + — , 

^ph / 



5p,p'5h,h' + -N 



relate the amplitude functions: 

Sc = C-6v . 

In the driving term on the l.h.s. of ( p.45| ) we use po,p/i = (C-p^)o,ph to obtain 

2 jd^r /iext(r; co) po,ph{r) = 2 C ■ 



,ext 



(3.51) 



(3.52) 



(3.53) 
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where the vector h*^"^* is built with the non- interacting states (c/. Pq^^^ in ( p.lOp ) 

Kphi^) = ^ph {h\ /ioxt(r; u) \p) . (3.54) 



Defining the cu— dependent matrices 

{hu+if] - eph)Sp^p:Sh,h' 



w 







the equations of motion assume supermatrix form [|| 



1 1 

2 2 ^ ' 



(5v = 2C ■ h' 



ext 



We now formally define a new, energy-dependent interaction matrix Vp_h(co') by 



o + -rjN + -N^^ - w 

2 2 



c 



c . 



Thus the response equations take the simple TDHF form 



n - Vp_h(a;) 



■5c = 2h' 



ext 



(3.55) 



(3.56) 



(3.57) 



(3.58) 



With this, we have reformulated the theory for a strongly interacting system in the TDHF 
form ( p.58| ) but with an energy dependent effective interaction. Our derivation has led to 
a clear definition of this effective particle-hole interaction and to a prescription on how to 
calculate this from the underlying bare Hamiltonian. 

The formal derivation appears to involve the calculation of the inverse of a huge matrix. 
The key point, however, is that the manipulation (|3.57| ) can be carried out diagrammatically. 
Then it becomes obvious that many terms occurring in the combination of matrices in (p.56|) 
are not part of Vp_h(i^)- Specifically, these are the chain diagrams in the direct channel [^. 



IV. DIAGRAMMATIC ANALYSIS AND LOCAL INTERACTIONS 



A. General strategy 



Generally, the non-local operators A/'(l, 2) and W(l, 2) in (p.l7 ) consists of up to 4-point 



functions. Cluster expansions and resummations have been carried out in Ref. 31 and led to 
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reasonably compact representations in terms of the compound- diagrammmatic quantities of 
the FHNC summation method. Nevertheless, due to their non-locality, it is difficult to deal 
with these quantities exactly. The simplest approximation for the operator is to keep just 
the local terms. These are given by the "direct-direct" correlation function Tdddri— r2|) of 
FHNC theory. This approximation is adequate but not optimal. 

On the other hand, summing A^o,pp'/i/i' over the hole states, Eq. ( [^.18| ), relates A/'(l,2) to 
the static structure function. Accurate results are available for S{q), either from simulations 



1 33, 34 1 or from the FHNC-EL summation technique [^8|, An alternative strategy to 
deal with non-local operators is therefore to demand that these results are reproduced in 
whatever approximate form one chooses to use. In this sense, by choosing A/'(l,2) to be 
local, naming the corresponding function rdd(r), and demanding that this operator in ( p.l8|) 
gives the known static structure function, we obtain the relationship 

%) = 5p(g)[l + fdd(g)5p(g)] (4.1) 

as a definition of rdd(f7) in terms of S{q). We adopt this view here and define the "best" 
local approximation for A/'(l,2) such that it reproduces the best known S{q). Since the 
exact S{q) contains a summation of exchange terms, this implies that their contribution to 
S{q) is mimicked by a local contribution to T^^^q). 

An "optimal" local approximation for the effective interaction W(l, 2) can be obtained 
along similar lines. From Eqs. (|2.14|) and ( p.ll|) we have 



H'o^pp'hh' — ^o,pp'hh' + - {^ph + ep'h') No^ppihh' ■ (4.2) 



The ground state Euler equation for pair correlations (|2.19|) implies that the Fermi sea 
average of -f^opp'h/i' vanishes. Postulating a local W(l, 2) ^ W{ri2), consistency relates this 
quantity to the local approximation of A^(l, 2). This leads to p8| 



mg) = -;|^fdd(g). (4.3) 



Our procedure of using the relationships ( |2.18D and (|2.19|) to construct local approxima- 



tions for No^pp/hh' and Wo^pp'hh' can be generalized to a systematic definition of optimal local 
approximations for the matrix elements of any non-local d— body operator: Averaging the 
matrix elements, which depend on d particle and d hole momenta, over the Fermi sea, gen- 
erates functions of the momentum transfers qj = p— hj only. Spelling out Fermi occupation 
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functions nh and np= 1 — np explicitly, this reads for a one-body quantity 
The extension to variables is obvious, 

d 



Cqi,-,qd - X] n jY^-' c)!\ Oo,pi...pahi...ha > (4-5) 

as is the extension to matrix elements Om,ny^o- 

We emphasize again that the quantities Oqj,...,q^ contain all exchange and correlation 
effects in a localized manner. Therefore, effects related to the Zph, as well as CBF corrections 
to the Cph , are already part of W{q) and TddiQ)- This implies, amongst others, 

Mp>h',ph ~ Sp^p'Sh,h' + {hp'lTddlph') , (4.6) 

and the relationship ( |3.51| ) between the supermatrices C and N simplifies to 

C = 1 + . (4.7) 

B. Matrix elements 

The localization procedure discussed above for M{1,2) implies 

N= ^fdd(g) ( '^'^'^ I npUpmhUi,, . (4.8) 

To simplify the notation, the 5q,±q' functions, together with the Fermi occupation numbers, 
are understood to be implicit in all the matrices from now on. Matrix products, i.e. sums 
over particle-hole labels, reduce to factors Sj,{q). The inverse of C is readily obtained from 
(13) as 

= 1 - 7^Xdd(g) I ^ ^ I • (4.9) 



2N 

with 



1 1 



X^q) = ^dd(g) _ 

l + 5p(g)rdd(g) 



In the spirit of the discussion in Sec. IV A , this is our definition of Xdd(Q')- According to 



( |A13|) , it can also be identified with the sum of all non-nodal diagrams. 
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Multiplying C ^ from both sides to ( p.57|) yields the u dependent effective interactions, 

V„_.m4(^^<'^ ^><'^ ^l). ,4.U) 

To summarize, the localization of A/'(l, 2) in an S{q) conserving manner has uniquely 
fixed the functions rdd(5') and Xdd((?) and, consequently, the corresponding matrices N and 
C^^. Calculating Vp„h(i^) from ( p.57| ) has thus been reduced to calculating V^^{q]uj) from 
W. 

In order to derive the explicit expressions, we need the optimal local form of ( p.43|) . This 
involves two steps, calculating the localized versions of the three-body vertices Kph^p'p"h'h" 
and K^p),f^,f^„ and deriving the inverse of the four-body energy matrix [E{ui)]~^. We 
expect these quantities to be sufficiently accurate within the convolution approximation, 
since improving on this only marginally changes the results [0 for bosons. 

The details of the derivation of the local three-body vertices Kq^qiqii and k'^^]^„ g defined 



in (|3.34| )-( |336| ) can be found in App. |B 3| . These are 



Kq^q'q" 



S{q')S{q") 



2m S,{q)S,{q')S,{q") 

s.w)s,{q'r - 

S{q')S{q") _ 
S{q')S{q' 



q- q'Xdd(g') + q- q"Xdd(g") - q'u3{q, q' , q" 



+ 



1 - 



(4.12) 



^rdd(g) 
4m 



SAq')SAq") 



q'Si'\q,q',q") 

SMSAq')SAq") 



+ 



q ■ q' ^ q ■ q" 



SAq') s^q") 



(4.13) 



Here, Sp {q,q',q") is the three-body static structure function of non-interacting fermions, 
defined in Eq. (P5|), and U3{q,q',q") is the ground-state triplet correlation function pB| . 
The implicit momentum conservation functions 5±q,q'+q/' ensure that both vertices depend 
on the magnitudes of the three arguments only. 

Going back to the Lagrangian, we realize that the term K^Jl„ r, is the coefficient function 
of the contributions to C'^^'^\t) containing Ui{t)U2{t) which we expect to be small. Our 
numerical applications to be discussed below will support this expectation. However, the 
vertex Kq^q'q" contains a term of the same form. Neglecting -R'^y/ q should, for consistency, 
also mean neglecting the same term in Kqq/gn which is then given by the very simple first 
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part of Eq. ( [4.12|) . In this term we recover, apart from Spiq) factors, also the Bose version 
of the three-body vertex. 



C. Effective interactions 



Next, the matrix elements ( |4.12|) and ([4.13|) are used in ( p.43| ) to calculate the dynamic 
parts of W, 



s. 



N 

q.-q' 
N 



W{q) + WA{q;oo) 
W{q) + WB{q;oj] 



(4.14) 



where the energy independent part W{q) has been defined in Eq. ([4.3| ). Because of the 
locality of the three-body matrix elements, we can write for the first dynamic contribution 

to (CT) , 



PiP2hih2 p\p'^h'^h'2 



Kp'^p'^h[hi2,p'h' 

piP2hih2,p{p2h[h'2 



K I I 

q,qiq2 ^^q^q2,q 



9191 



1^ E i^M-' 

hih2h'^h'2 



J piP2hih2,p[p'2h[h'2 



-^^^q^q^q2^ \qi,q2;uj)K, 



9192,9 



(4.15) 



9192 



with implicit factors 5q,qi+q2 ^q,q'^+q^ for momentum conservation. The other contributions 
to ( |3.43| ) are calculated analogously. The inverse four body energy matrix and the pair 
propagator 

(4.16) 



hh'h"h"' 



pp'hh',p"p"'h"h"' 

are calculated and discussed in App. |C|. Basically, the pair spectrum is built from two 
particle-hole spectra. These are, however, not centered around free particle spectra but 
around the Feynman dispersion relation. Consequently, our pair propagator also includes 
two-phonon intermediate states. 

The resulting expressions for the energy-dependent ^^^(g;^) are then 



WAiq;uj) 
WB{q;uj) 



1 



OAT Z-^ V ' 



2N 

1 

2iV 



'2 E-\q\ q"; u) + \KlJ>J' E-'\q\ q"- -co) , (4.17) 



'(9) |2 fn-l 



E f^?i",oW (E'\q\q';u;)+ E-'*{q\q"; 



(4.18) 



q'q" 
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Similar to the boson theory, the dynamic parts of the interactions are expressed in terms 
of three-body vertices and an energy denominator, the latter now being "spread" over the 
whole width of a two-particle-two-hole band. 

The last step in our formal derivations is the calculation of Vp_h('^)- Carrying out the 
operations ( ^.57] ) yields the energy- dependent, but local functions 

K(g;u;) = \/p_h(g) + [<p u;) + [a-]^ ^^(g; -a;) 

+ {W,iq;uj) + W:iq;-U)) , (4.19) 

K(g;a;) = Vh(?) + [a^f W^iq; u) + [a;]' W* (q; -u) 

+ (W,iq;u^) + W:iq;-Uj)) , (4.20) 

with a± ^ [S,{q) ± S{q)]/2S{q) . 



V. DENSITY-DENSITY RESPONSE FUNCTION 



A. General form 



We now derive the density-density response function x(g;a;). The final result for the 
dynamic effective interactions, (|4.19|) , ( [4.20| ), is inserted into (|3.58| ), which is solved for 5c. 
The induced density is then obtained from Eq. ( p.47| ). Using Po,p/i('^) ~ j^e'^^^'^^^^^ we obtain 



h 

^^^^^[6c^'^\q;u;) + 6c^'-\q;oj)] , 



^Kphi^) = (±(^^ + i^) -eph)^cg^(w) 

- V {q-uj)5c^^\q-uj) - V:{q--uj)5c^^\q-u) 



(5.1) 



where we abbreviate in the second line 6c'^^^\q]ijj) = -^^h^'^pit\^)- Spelling out 
Eqs. dpI) explicitly. 



(5.2) 



dividing by {±{JTjjj + iri) — Cph) and summing over h yields 



6c^^^\q; 



CO 



I KAq- u) + V^{q; co) 6c^'^\q; co) + V;{q; -u) 6c^'^\q; u) 



Ko(g; w) 



(5.3) 



23 



with the positive-energy Lindhard function 

1 

N 



(5.4) 



(5.5) 



^ hw -eph + iri 
which is related to the full Lindhard function by 

Xo{q; w) = Ko(g; uj) + kKq; -uj) . 

Solving for 5c''^'^\q; to) and inserting into we obtain for x(g; u) 

X{q;uj) = N{q;uj)/D{q;uj) 
N{q;uj) = KQ{q;uj) + KQ{q;-uj) 

D{q;uj) = 1 - KQ{q-,uj)V^{q-,Uj) - Kl{q-,-uj)V*{q-,-u) 

+ KQ{q;uj)K*o{q;-uj) Vj^{q] u)V* {q; -u) - V^{q-, Uj)V^* {q; -uj) . (5.6) 

Eq. (|5.6| ) is the TDHF response function for local and energy dependent interactions. Ev- 
idently, the conventional RPA form ( |1.7| ) can only be recovered if the interactions V^{q]u) 
and Vg{q;u) are energy independent and equal. Clearly, our result ( |5.(j| ) significantly dif- 
fers from ( |1.7| ) with V^-ilq) simply replaced by some energy dependent Vp_h(g;a;). Such an 
RPA- like form for the density-density response function lacks microscopic justification. 



B. Long wavelength limit 

In the limit g — )■ 0, the spectrum is dominated by collective excitations, e.g. zero sound or 
plasmons. Both vertices ( [4.12| ) and ( |4.13|) vanish linearly in q, hence W^(g;a;) and W-^{q;u) 



are quadratic in g as g — 0. 

For neutral systems, the dynamic corrections to the effective interactions V^^{q;u) in 
( 4.19| ), ( 4.20 ) are therefore negligible in the long wavelength limit. The long wavelengths 



density-density response function is then given by its RPA form ( |1.7] ), with the static particle- 
hole interaction Vp-h(g). The zero sound speed Cq is determined by the long wavelength 
solution of the RPA equation. 

For charged quantum fluids, cr^ ~ Sp{q)/2S{q), hence V^{q,u) = V^{q,u), which again 
implies the RPA form ( |1.7| ) 

X{q;(^) = -, -F7-, r as q^O. (5.7) 

^- Xo{q](^)Vj,{q;uj) 
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However, now the effective interaction is 



V^{q;uj) = V^4q) + 



W, {q;io) + W, (g; -oo) + W^{q;io) + {q; -u) 



AS\q) 

The static particle-hole interaction approaches the Coulomb potential Vc{q) = 47re^/g^ 



as g — )• . 
(5.8) 



Vp-hiq) = Vciq) + Vo as g^O . 



(5.9) 



We can therefore write (B78|) as 



Va{Q',^) = V^{q;uj) = fc(g) + Vo(w) as q-^0. 



(5.10) 

As for charged bosons ||36|, the two-pair fluctuations modify the RPA result. The static 
potential Vp-J^q) and Wa,b{(1'i<^) contribute for g— t-O at the same level. 



C. Static response function 



E ^(g, g'; a; = 0) is real and negative, this is most easily seen from the representation (|C9|). 
Therefore, all interactions W^^{q;0) in (^T?D-(^:i8l) and V^,s{q;0) in (|CT|)-(|i:20|) are real. 
The response function ( [5.6| ) can again be cast into the RPA form 

xo(g;0) 



with a static effective interaction 

^tat(g) = K>-h(g) + 



i-KtM Xo{q;0) 



2S'{q) 



W;(g;0) + W;(g;0) 



(5.11) 



(5.12) 



Unlike Eq. ( ^.8|) , this form holds for all wavelengths. 



For short wavelengths the static response function has the asymptotic form B^, |38 



x(g-)-oo;0) 



(T) 



+ 0{q- 



(5.13) 



t{q) 3t2(g) N 

where (T) is the kinetic energy. In the RPA, one obtains in Eq. ( |5.13|) only the kinetic 
energy of the non-interacting system. To obtain the correct asymptotic form, it is therefore 
necessary to include pair and, possibly, higher order fluctuations. 

Again, we know the result for bosons as a guide: treating pair fluctuations in the "con- 
volution" approximation leads to the correct asymptotic behavior with (T) in (|5.13 ) given 



in that approximation |18 
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We show in App. that 



1 2 (T) T 

Ktat(g->oo) = -W;(g^oo;0) = , (5.14) 



where (T)'-'^ is the kinetic energy in "uniform hmit" or "convolution" approximation ( |A12 ) 



Hence, inserting the short wavelength expansion of the Lindhard function, the static response 
function (|5.11|) indeed assumes the form (|5.13|) 



\CA 



/ , 2 8 (Tf 

X(g;0) = -^-^^ as g->oo, (5.15) 



with the kinetic energy being calculated in the uniform limit approximation (|A12|) 



VI. APPLICATIONS 

A. Dynamic structure of '^He 

1. Motivation 

The helium fluids are the prime examples of strongly correlated quantum many-body 
systems. They have been studied for decades, and still offer surprises leading to new insight. 
It is fair to say that understanding the helium fluids lies at the core of understanding other 
strongly correlated systems. The most important and most interesting field of application 
of our theory is therefore liquid ^He. 

Recent developments ^ have brought manifestly microscopic theories of ^He to a level 
where quantitative predictions of the excitation spectrum are possible far beyond the roton 
minimum without any information other than the underlying microscopic Hamiltonian ( p.ll ). 
^He is the more challenging substance for both, theoretical and experimental investigations. 
Experimentally, the dynamic structure function S{q; u) of ^He is mostly determined by 
neutron scattering. The results are well documented in a book the theoretical and 



experimental understanding a decade ago has been summarized in Ref. ^ Recent inelastic 
X-ray scattering experiments have led to a controversy on the evolution of the zero sound 
mode at intermediate wave- vectors |]4ll-43|, we will comment on this issue below. 



The RPA (13) suggests that S{q; u) can be characterized as a superposition of a collective 
mode similar to the phonon-maxon-roton in ^He, plus an incoherent particle-hole band 



which strongly dampens this mode The picture is qualitatively adequate but misses 
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some important quantitative physics: In ^He the RPA, when defined through the form ( [L.7|) 
and such that the sum rules ([L.4[)-(pr5D are satisfied, predicts a zero-sound mode that is 
significantly too high. This is consistent with the same deficiency of the Feynman spectrum 
( |1.6D in ^He. Drawing on the analogy to ^He |Q, the cure for the problem is, as pointed 
out above, to include pair fiuctuations ^u^^,f^y{t) in the excitation operator. 

An alternative, namely to lower the collective mode's energy by introduction of an effec- 
tive mass in the Lindhard function, leads to various difficulties: First, one violates the sum 



rules (|1.4| )- (|1.5| ), i.e. one disregards well established information on the system. Second, 
the effective mass is far from constant; it has a strong peak around the Fermi momentum 



45-48 , a secondary maximum around 2A;f5 and then quickly falls off to the value of the 



bare mass. In fact, it is not even clear if the notion of a "single (quasi-)particle spectrum" 
that is characterized by a momentum is adequate at these wave numbers. 

The localization procedure of Sec. |^ implies that the only input needed for the applica- 
tion of our theory is the static structure function S'(g), whereas the single-particle spectrum 
is that of a free particle. We hasten to state that we do not claim that the precise location 
of the single-particle spectrum is completely irrelevant for the energetics of the zero sound; 
we only claim that the dominant mechanism in Bose and Fermi fluids is the same, namely 
pair-fluctuations. In order to maintain the sum rules ( |1.4| )-( pr^ ), any modiflcation of the 
particle-hole spectrum must go along with an inclusion of exchange effects. At the level of 
single-particle fluctuations 0, ^, such a calculation is quite feasible [|^, However, to 
describe the dynamics of ^He correctly, it is insufficient to include only the CBF single par- 
ticle energies ( p. 141) . These suggest a smooth spectrum with an effective mass slightly less 
than the bare mass, in contradiction to the highly structured spectrum mentioned already 
above. 



2. Collective mode 



For our calculations we have used input from the FHNC-EL calculations of Ref. ESI that 



utilizes the Aziz-II potential and includes optimized triplet correlations as well as four- 
and five-body elementary diagrams. An overview of our results for bulk ^He and a compari- 
son with both the RPA and experimental data is shown in Figs. |I] for four different densities. 
The most prominent consequence of pair fluctuations is a change in energy and strength of 
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FIG. 1: (Color online) S{q;uj) of ^He, for the densities p = 0.0148, 0.0166, 0.018, 0.021^3. The 
experimental results for the collective mode (dots) are from inelastic neutron scattering experi- 
ments at the ILL (Ref. ||). The densities 0.0166, 0.0180 and 0.0200 correspond in good 
approximation to the pressures p = 0, 5, 10 bar ||5^, |5^]. Dashed lines are equidistant contours 
marking the same absolute value in all plots. Solid lines are the boundaries of the particle-hole 
continuum for m* = m. The blue boxes show the RPA result for the collective mode. 

the collective mode and its continuation into the particle-hole band. Pair fluctuations also 
contribute a continuum background outside the particle-hole continuum. 

At long wavelengths, the collective mode is sharp and well defined above the particle-hole 



band, exhausting most of the sum rules (1^) and (]_^). In this regime, the RPA provides 



a faithful description of the physics. This is in accordance with the observation that the 
dynamic correction to the effective interactions vanish, for neutral systems, in the long- 
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wavelength limit. With increasing density, the speed of sound increases and the phonon 
becomes farther separated from the particle-hole band. 

Further details are shown in Fig. 0. At intermediate wavelengths the collective mode 
bends down due to the attractiveness of the effective interaction. This is where the dynamic 
theory starts to deviate visibly from the RPA. Evidently, pair fluctuations are the major 
cause for lowering the energy of the collective mode, although they do not completely bridge 



the discrepancy between the RPA and experiments |2J, |2^. This is expected because, for 
bosons, pair fluctuations bridge only about two thirds of the gap between the Feynman and 
the experimental roton energy |17|. Three-body and higher-order fluctuations are also 
important [18|. We expect that these corrections are smaller in ^He due to its lower density, 
yet not negligible. 

When the collective mode enters the particle-hole band, a slight kink in the position of the 
maximum in S{q; u) is expected, as well as an abrupt broadening of the mode. At saturated 
vacuum pressure, shown in the left part of Fig. |^, these effects are difficult to identify in the 
experiments |2^. A possible reason is that the observed mode stays always very close to 
the particle-hole band. The measured mode width in Fig. |^ gives no clear indication of the 
upper boundary of the particle-hole band other than that it is determined by a spectrum 
with an average effective mass of m* ^ m. 

The situation is much clearer at higher pressure: With increasing density, the speed of 
sound increases, separating the collective mode farther from the particle-hole band. For 
p= 0.02 A~^ a clear kink is identified at g~ 5nm~^ (Fig- il right part). The broadening is 
also more abrupt and, in particular, does not increase for larger values of q. Similar to SVP, 
explaining these data requires a boundary of the particle-hole band that is even above that 
of the non-interacting Fermi fluid. Damping due to multiparticle excitations is, on the other 
hand, for both densities far too small to account for the experimentally seen broadening of 
the zero sound mode. 



3. Frequency dependence of S{q;uj) 

For a quantitative discussion we show in Fig. ^ the dynamic structure factor as a function 
of frequency at a sequence of wave vectors. We conclude that the RPA quantitatively and 
even qualitatively differs from our theory and the experiment. Including pair fluctuations 
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FIG. 2: Zero sound mode calculated within the pair fluctuation theory (full blue line), RPA (red 
chained line) and experimental data by the ILL group [^] (square symbols) and (circles). The 
bars indicate the width of the fit to the data, the line at the bottom of the figure gives the width 
due to pair fluctuations enhanced by a factor of 10 to make it visible. The dashed blue line gives 
the FWHM of the mode within the particle hole continuum. Left part: p = 0.0166A~'^, right part: 
p = O.O2A-3. 

improves the agreement with experiment significantly. The arrows in panes (c) and (d) 
indicate the maximum of the experimentally observed dynamic structure function. 

In Fig. ^(b) we also show the consequence of the plausible simplification of our theory 
discussed already in connection with Eqs. ( ^4.12| ) and ( 4.131 ): We neglect all terms that 



vanish for bosons as well as for large momentum transfers q,q',q" > 2kp. This is K^jln q 



and, consequently, the second term in Kqq/q//, Eq. ( 4.12 ). The three-body vertex is then 
given by the first term in Eq. (|4.12|) , see also (p 1|) . This simplifies the effective interactions 
significantly: Only the first term of Eq. ( [4. 171) for WA{q]u) contributes, and Wsiq^u)) is 
neglected. Fig. |^(b) shows that these simplifications modify our results only marginally, the 
form (|D1|) can therefore be considered a practical and useful simplification of our theory. 

Figs, ^(c) and |](d) show our results for the two momentum transfers q = 2.4 kp = 
1.89 A^^ and q = 3.2 = 2.52 A^^. Recent X-ray scattering experiments in that momentum 
range [^T|-^ appeared to support the notion of a high-momentum collective mode without 
visible damping by incoherent particle-hole excitations. Figs. |^(c) and ^(d) show that pair 
fluctuations lead to a narrowing of the strength of S{q] u) compared to the RPA. To facilitate 
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the comparison with experiments, we have convoluted our result with the instrumental 
resolution of 1.58 meV, the results are also shown in Figs, ^(c) and 0(d). After this, our 
results agree quite well with the experimental spectrum. Also, the location of the observed 
peak intensity for q = 2.4 kp appears to be consistent with our calculation. The RPA is, 
on the other hand, too broad to explain the data. We also point out that a value of the 
effective mass close to m* ^ m is consistent with our theoretical calculations . We have 
to conclude therefore that the observed width of the X-ray data are also consistent with our 
picture. 

After a regime of strong damping we see in Figs. an intensity peak at momentum 
transfer of g ~ 2.5A;f- With increasing density, this peak moves towards the lower edge 
of the particle-hole band and becomes sharper. Such a peak should be identified with the 
remnant of the roton excitation in "^He, broadened by the particle-hole continuum. The 



overall agreement with the experiment is quite good, see Fig. 1 of Ref. |2^. Our theory 
predicts a "roton minimum" that is slightly above the observed energy; this is expected 
because for bosons a similar effect is observed. To obtain a higher accuracy, triplet- and 



higher order fluctuations must be included |18 



4- Static response 

For completeness, and because the quantity should be obtainable by experiments and 
simulations similar to those for ^He [^4|, ^ and on bulk jellium we show in Fig. ^ the 



static response function 0) of ^He at p = O.OIGGA-^. The main peak, which is a result 
of the local symmetry in the fluid, is visibly raised compared to the RPA result. We suspect, 
form experience with the boson theory, that this peak is still a bit underestimated. 

The comparison also lets us assess the validity of an energy independent particle hole 
interaction. Fig. ^ shows a comparison between the FHNC Vp-Y,{q) and the static effec- 
tive interaction ( [5.12| ). Evidently, the qualitative structure is very similar, in particular 
K-h(Q'~^0) — ^statil — > 0) as discussed in Sec. |V B| . The most visible difference is that 
VstntiQ) approaches a constant for large q, see Eqs. (|5.13|) and ( |5.14D . 
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FIG. 3: (Color online) S(ci\uj) for ^He as a function of energy at /) = 0.0166A^^ for a sequence 
of momentum transfers q = 0.8, 1.6, 2.4, 3.2 kp (a)-(d). Also shown is the RPA (dashed, red). 
The solid blue line is the result of this work with the simplified Wj^{q;uj) and W^{q,uj) = as 
discussed in the text. In pane (b), we also show the results when the full Wj^{q;uj) and W^{q;uj) 
of Eqs. ( 4.17| ) and ( [4.18 ) are retained (short dashed magenta line). The results from the different 
approximations are almost indistinguishable in panes (a),(c) and (d) and therefore not shown. The 
black dash-dotted line in panes (a) and (b) are fits to the experimental results of Ref. In panes 
(c) and (d) we indicate the maximum of the experimentally observed dynamic structure function 
by an arrow. We also plot in panes (c) and (d) recent inelastic X-ray diffraction data obtained 



by Albergamo et al. |41] (boxes) as well as our theoretical results folded with the experimental 
resolution (dashed line). 
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FIG. 4: (color online) Static response of ^He at p = 0.0166A-3. The red curve shows the RPA 
result whereas the blue line is the result of this work. 
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FIG. 5: (color online) Effective interaction of ^He at /? = 0.0166A ^. The red curve shows the 
static effective interaction Vp-i,{q) whereas the blue line is Kitat('7)- 

B. Electron liquid 

The second typical area of application of microscopic many-body methods is the electron 
liquid It provides the basic understanding of valence electron correlations in simple 



metals. In its two-component version it has proved useful for describing the electron-hole 
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liquid in semiconductors. 

Compared to the helium fluids, the soft repulsion of the Coulomb interaction induces 
substantially weaker correlations. Therefore, electrons are much less challenging than ^He 
and the RPA (or slightly modifled versions) contain much of the relevant physics. 

Correlations are somewhat more pronounced in layered realizations of the electron liquid, 
such as Si- and GaAs-AlGaAs hetero-structures. For electrons on He surfaces preliminary 



results show that at very low densities, again, a roton-like structure evolves for inter- 
mediate wave vectors. 

We have seen that pair fluctuations contribute, already at long wave lengths, to the 
static response function, see our discussion in Sees. [VB| - |V Q . Most important are, of course, 
those effects that are qualitatively new consequences of multiparticle fluctuations. These 
are the short-wavelength behavior of the static response function and the appearance of a 
new feature in the dynamics structure function, namely the "double plasmon" excitation. 



The latter has raised new interest ||26| , p7| in studying the dynamics of electrons at metallic 
densities in this (g; u) region. 



1. Double Plasmon 

Figure |^ shows the dynamic structure factor S{q; u) obtained from the pair fluctuation 
theory. We have chosen two different densities p = 3/(47rr^a|), corresponding to Al, = 
2.06, and Na, = 3.99. Immediately obvious are the flnite width {i.e. lifetime) of the 
plasmon above the particle-hole band, and a second peak-like structure around twice the 
plasma frequency Up. 

Characteristic cuts at constant wave vectors q are shown in Fig. ^ for Na. In parts (a) 
and (b) the plasmon is outside the particle-hole band and rather sharp; the second peak 
slightly above 2hup = 4.5 tp is clearly visible. We identify this feature, which has also been 
observed experimentally [^, with the "double- plasmon" . 

The "double-plasmon" excitation is due to the emergence of an imaginary part in V^{a, u) 
at u = 2up, caused by the appearance of an imaginary part of the pair propagator 
{q' , q" ] u) . It is therefore a genuine multipair effect. The properties of the pair propa- 
gator are discussed in in App. |C 2| . From (|C19| ) we obtain for the double-pole part of the 
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FIG. 6: (Color online) The figure shows S{q;uj) of an electron liquid with density parameters 
= 2.06 and = 3.99 appropriate for Al and Na, respectively. As in Figs. |l], dark red regions 
correspond to high intensity (logarithmic scale). The blue line is the position of the double- 
plasmon peak obtained in the present work, red dots are experimental results from inelastic 



X-ray scattering and green diamonds results from Green's functions calculations [27 
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In Fig. 0(c), the plasmon is broad and Landau-damped, while the double-plasmon still 
shows a clear structure, even at the brink of entering the particle-hole continuum. Some 
structure in the spectrum persists to even higher momentum transfers: At g = 2.0 k-p in 
Fig. 0(d), traces of the ordinary as well as the double plasmon show up as a faint double- 
peak structure, with its minimum where the RPA yields a single maximum. 

We now investigate the nature of the slight but measurable peak in the loss function 
at approximately twice the plasmon frequency Wp. Fig. || shows S{q,u) for = 3.99 for 
three different momentum transfers, the position of the double plasmon is marked with 
arrows. 

We have already shown in Figs. |] the location of the double plasmon excitation and 
a comparison with the experimental inelastic X-ray scattering data ^6, 27 . The double- 



plasmon is also accessible by Green's function methods These results are very close to 
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FIG. 7: (Color online) S{qo;uj) for Na (r^ = 3.99), at wave vectors qo (a) 0.15 fep, (b) 0.6 A;f, (c) 
1.3 kp, and (d) 2.0 kp. The full (blue) lines are our pair fluctuation theory, dashed (red) lines are 
the RPA results using Vp-h{q)- To make the plasmon visible, the RPA data have been broadened 
artificially by adding an imaginary frequency of 10~^eV/?i. The dotted (green) lines in (a) and (b) 
refer to neglecting K^Jgii q in Eqs. ( [4.12| )-( |4?T^ ), and the dash-dotted (black) lines include ground 
state triplet correlations. At larger momentum transfers these effects are too small to be visible. 



those of our pair fluctuation theory. This can be understood from the fact that the leading 
terms of the long-wavelength part of the pair propagator actually contain no correlation 
effects, see Eq. |C29| . Hence, theories that are less well suited than CBF for the description 
of strong correlations should, similar to the single plasmon, give the right answer. The 
remaining discrepancy with experiments must therefore be attributed to lattice effects. Fig. 
^ shows more details of S{q,u) at a sequence of three different momentum transfers for 
Ts = 3.99 (the position of the double-plasmon is marked with arrows), in particular in order 
to assess the relative strength of the double-plasmon excitation compared to the underlying 
continuum. 
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FIG. 8: Cuts of the density-density response function at Na-density (r<j = 3.99), for constant 

momentum transfer q = 0.4 /cp (solid line), q = 0.6 kp (dashed line) and q = OAkp (dash-dotted 
line). The arrows mark the position of the double-plasmon. 

2. Static Response 




FIG. 9: Static response function (left), and static effective interaction (right) of the electron liquid 
at r<j = 3.99. Full blue lines are our results, black dash-dotted lines a fit based on the simulations 



56, pfl]. Dotted red and thin broken lines show the RPA with Vp-h{q) and Vc{q), respectively. 



Monte Carlo studies of the static response function x(g; 0) were performed for two- and 
three-dimensional ^He [^4|, ^ and on bulk jellium for = 2, 5 and 10. While x(g; u) 
is accessible experimentally, for electron liquids it is popular to define a static local field 
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correction to the Coulomb interaction Vc{q) via 



^tat(g) = d,{q) {1 - G{q)) . (6.2) 

From our analysis it is clear that a response function in the RPA form can be defined only 
for g — > and at u = 0. Therefore, only in these two cases such a function is a physically 
meaningful quantity. 

In the q ^ CO limit, our theory yields a finite value for Ktatiq), resulting in G{q) oc g^, 
whereas Vp_h('?) falls off like the bare potential. This correct g— dependence arises solely from 
multiparticle fluctuations. In Fig. ^ we compare our results with the Monte Carlo data, and 
with curves calculated from an analytic analytic fit for —Vc{q) G{q) obtained from the latter 
||60| . The agreement is remarkably good. 

No trace of a possible "hump" in G{q) around 2kp as a remnant of some charge- or spin- 
density wave instability was found in the simulations, but it also was not fully conclusively 
ruled out. Our results, clearly, do not yield any such peak structure at 2A;f either. 

VII. SUMMARY 

We have presented the fermion version of theories of the dynamic response of Bose flu- 
ids that have been developed in the past successfully by Jackson, Feenberg, and Campbell. 
These methods form the basis of our present understanding of the dynamics of Bose fluids. 
Our derivations were admittedly lengthy but eventually led to a reasonably compact formula- 
tion of the dynamic response of correlated Fermi fluids. Our final result could be formulated 
as a set of TDHF equations in terms of dynamic and non-local effective interactions. 

For the first applications we have reduced the theory to a practical level capturing the 
relevant physics, while avoiding many of the technical complications. In particular the 
version of the equations of motion spelled out in Appendix ^ has proved to be adequate for 
systems as different as ^He and homogeneous electrons. It is hardly more complicated than 
TDHF. The sole required input is the static structure function S{q) which can, in principle, 
also be obtained from simulations. Our developments have led to quantitative improvements 
of our understanding of ^He and electrons as well as to the description of qualitatively new 
effects like mode-mode coupling, multiparticle spectra, and damping. 

We have, at various places, commented on the role of the particle-hole spectrum. In 
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the homogeneous electron hquid, the interaction corrections to the single-particle spectrum 



are relatively small p5| , |6T| , the theory formulated here should therefore suffice for many 
purposes. The situation is more difficult in '^He: As is seen from our results, good agreement 
with experiments can be reached by assuming a spectrum of non-interacting fermions. In 
particular looking at the zero-sound damping suggests that, at q ~ kp, the boundary of 
the single-particle continuum should be close (perhaps even above) to the one given by a 
non-interacting spectrum, cf. Fig. ^. This is not in contradiction to experiments [p2i |62| 
suggesting an effective mass ratio m*/m 3 at the Fermi surface. One reason is that 
the effective mass ratio drops rapidly with distance from the Fermi surface. The more 
fundamental reason however, is that the concept of describing the particle-hole excitations 
by a spectrum that depends on momentum only is questionable at elevated wave numbers. 
More precisely, the single-particle motion is described by a non-local, energy dependent self- 
energy. Upon closer examination it becomes clear that exchange effects are intimately related 
to self-energy corrections and exchange effects must therefore be included simultaneously. 

In independent work, we have used the ideas of CBF theory as well as the Aldrich- 
Pines pseudopotential theory to calculate the single-particle propagator in '^He. In both 
three and two dimensions, we found good agreement between the theoretical effective mass 
near the Fermi surface, and that obtained experimentally from specific heat measurements 
J^, ^ However, the somewhat ad-hoc use of the effective interactions in that work is 



still awaiting rigorous justification. This is the subject of future work. 
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Appendix A: Ground state theory 



The essence of FHNC-EL 



For the sake of the discussions of this work we here briefly review the essence of variational 
FHNC theory. The diagram expansion and summation procedure that is used to derive, for 
the variational wave function ( p.2| ) a set of equation for the calculation and optimization 
of physical observables has been described at length in review articles [2T| and pedagogical 
literature [^. Details on the specific implementation for ^He are given in Ref. 

Here, we spell out a reduced set of equations. These do not provide the quantitatively best 
implementation p8|] of the FHNC-EL theory, but they contain the relevant physics: They 
provide, in the language of perturbation theory, a self- consistent approximate summation 
of ring- and ladder diagrams ||2^, thereby capturing both, long- as well as short-ranged 
features. 

In the simplest approximation which contains, as we shall see momentarily, the 
"RPA" expression ( |1.7|) , the Euler equation ( |2.5| ) can be written in the form p8[ 

SF{q) 



Siq) 



(Al) 



where t{q) = h q /2m is the kinetic energy of a free particle, and 



V,4r) = [1 + rdd(r)] vir) + - V^/T+f^ 



m 



+ rdd(r)wi(r) 



(A2) 



is what we call the "particle-hole interaction". Auxiliary quantities are the "induced inter- 
action" 



wi{q) = -t{q) 



n 2 r 



Spiq) S{q) 
and the "direct-direct correlation function" 



S{q) 
Spiq) 



(A3) 



fdd(g) = {S{k) - SM)/sUq) 



(A4) 



(see also Eq. ( [4.ip ). Eqs. ( [A1[ )-([A^) form a closed set which can be solved by iteration. 
Note that the Jastrow correlation function has been eliminated entirely. 

The relationship ( [Al| ) between the static structure function S{q) and the particle-hole 
interaction Vp-J^q) can also be derived from Eq. ( p. . 7p , if the Lindhard function is replaced 
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with its "mean spherical" or "collective" approximation (CA), 

2t{q) 



CO 



(A5) 



(hco + ir]y-t^{q)/SUq) ' 
The essence of this approximation is to replace the branch cut in xoiQi^) by a single pole; 
its strength chosen such that the first two sum rules agree when evaluated with the full 
Lindhard function xoil^cu) or in the collective approximation Xo^iQ',^), i-^- 



I duj Xq^Wi^) = duxoil'-,^) 



du u Xo^ic[',co) = 3m du u Xoil] cu) 



(A6) 



In fact, ( |1.7| ) together with (|A5| ) or, alternatively, 

t{q) ( 1 



2 \S\q) Sl{q) 



(A7) 



can be used |28| to define the particle- hole interaction from an accurately known S'(g). 
The energy, consisting of kinetic and potential energy (T) + (V), is 



E = -Ntp + Er + Eq 





Er 
Eq 



N r d^q 



m 



t{q) [SUq)~l-S\q) + S{q)] T^q) ■ 



4 J (27r)2p 

Here, tp is the Fermi energy, and, in this approximation, 

c{q) = sM-i + {sUq)-i)rM. 



(A8) 
(A9) 

(AlO) 
(All) 



To make the connection with the limiting behavior of xili 0) in Sec. |V (J| , we next spell 
out what is known as the "uniform limit" or "collective" approximation (CA). Products of 
functions which in coordinate space vanish for r — )■ oo are considered small. This implies to 
expand Va/T+T^mIO ~ |Vrdd(^') and to neglect C(r). The kinetic energy then is 

= Tf + \Y1 ^(^) ^(^) • (A12) 

q 

Here, Tp = 3Ntp/5, and Xdd(q') is the "non-nodal" function. In our reduced FHNC approx- 
imation, Xdd{q) is related to the static structure factor by 

1 1 



Xdd(g) 



Sp{q) S{q) 



(A13) 



41 



P h 

V 



P h 

V 



o 



P h 

V 
4 




6 



V 



FIG. 10: Diagrammatic representation of some contributions to po,p/i(i")- The upper row shows 
the diagrams defining the local approximation. The second row are the leading exchange diagrams 
and the third row shows two corrections due to the non- locality of 7\A(1,2). 

Appendix B: Diagrammatic analysis 

1. Transition density 

We first examine the diagrammatic structure of CBF matrix elements Pophi^) of the 
density operator, ( |3.10| , pm| ). The simplest approximation for Mph^pih' has been spelled out 
in Eq. (|4.(j|), the corresponding approximation for Po,p/i(^) 



Po,p.(r) = pipdr)+P d'r' d'r" 6{r-r' 



r-r'|fcp) rdd(r'-r")p^,,.(r"). (Bl) 



The diagrammatic representation of some leading diagrams contributing to Pophir) is 
shown in Fig. 0. As usual, open points represent particle coordinates rj, while filled points 
indicate an integration over the associate coordinate space and a density factor. Dashed 
lines connecting points and tj represent a function Tddifij), and oriented solid lines an 
exchange function i{rijkF). New elements are particle- and hole-states, depicted as upward 
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FIG. 11: Diagrammatic representation of the local approximation for Mph^p/hi. 

(particles) or downward (holes) lines entering or leaving the diagram. 

The three leading terms ( |B1| ) are shown in the upper row of Fig. TU. In the second 



row of Fig. ^ we show the leading exchange diagrams. In the representation (|3.1CI|) , these 



originate from the factors Zph in the definition of the po^ph{r), these are shown as the first 
two diagrams. Exchange terms also originate from the matrix element {ph'\rdd\hp')^, these 
are shown as third and fourth diagram in that row. Evidently there is a partial cancellation. 
The diagrams shown in that row also serve as an example for how the representations (p.lO|) 



and (|3.11|) are equal: Starting from the form (|3.11| ), the diagrams originating from the Zph- 
factors {i.e. the first two diagrams in the second row), have opposite signs; and the exchange 
term of {pp'\^dd\hh')^ yields the third diagram with interchanged particle- and hole labels. 
The sum of all three diagrams is the same. 



2. The M(I) matrix 

Our next task is to show that the diagrams representing M'^l^i^ni^if^,, are a proper subset 
of those contributing to Mph^p>p"h'h"- We restrict ourselves here to the simplest case, which 
is the numerically implemented version. We start with the two-body matrix Mph^p>h' ■ As 
spelled out in Eq. (|4.6| ), besides the 5-function, the leading contribution is the local term in 
the two-body operator 

Moe(l,2) = rdd(ri2). (B2) 

The diagrammatic representation of this approximation for Mph^p>h' is shown in Fig. |Tl|. 

A diagrammatic expansion of the matrix elements Mph^p'p"h'h" can be derived in exactly 
the same way as the corresponding expansions of the two-body matrix elements [^. Gen- 
erally, the Mph^p'p"h'h" are matrix elements of a non-local three-body operator, which can 
be expressed in terms of FHNC diagrams. Restricting ourselves again to the numerically 
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implemented level, we need these quantities in an approximation equivalent to the "uniform 



limit approximation" [|17| for bosons. We generalize this approach to fermions by keeping all 
diagrams contained in the Bose case plus those, where the end points of the correlation func- 
tions are linked by exchange paths (the bosonic g{rij) —1 is identified with the direct-direct 
correlation function Tdd{rij) ). This procedure has already been used for deriving the optimal 
triplet correlations for the fermion ground state The diagrammatic representation of 



this approximation is shown in Fig. ^ the analytic form is 

= Kh' (p/^"|rdd(i,2)bV') - V (/i'/i"|rdd(i,2)|V') 
1 

'2 



+^(p/i'/i"|rdd(3,i)rdd(i,2)|V/) 



hi hi 

+ (p/i'/i"|rdd(i,2)rdd(2,3)|V/) 

- J2{ph'\^M\hhi) {h"h\Tdd\p"p') - 5^(p/ii|rdd| V) {h'h"\TdA\hp") 

hi hi 

+ {ph'h"\T^,\{l,2,3)\hp'p") 

+ {{P'h') ^ {p"h")} . (B3) 

Here, in convolution approximation, 

rdd1i(ri, r2, rs) = ^ j d^riVddiri - r4)rdd(r2 - r4)rdd(r3 - r4) 

+ — d^Vi d^r^ £^(|r4 - r5|fcF)rdd(ri - r4)rdd(r2 - r5)rdd(r3 - rs) 
2z/ J 

^ ~ j d^rid^r^l^{\ri - r5|A;F)rdd(ri - r4)rdd(r3 - r4)rdd(r2 - rs) 

+ ^ y"c/V4rfVsc?V6£(|r4-r5|fcF)^(|rs-r6|A;F)^(|r6-r4|A;F) 

xrdd(ri - r4)rdd(r2 - r5)rdd(r3 - re) . (B4) 

The first two lines are invariant under exchanging r2 -H- rs, equivalent to exchanging {p'h') o 
{p"h") in (§1). 

Optimized triplet correlations improve the description of the ground-state structure, in 
particular in the area of the peak of the static structure function and also improve, for 



bosons, the density dependence of the spectrum ||T^. These correlations add another term 



to the three-body function rj^(ri, r2, r^). The expressions are lengthy |2^, we refrain from 
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spelling them out here and just show the diagrammatic representation of some typical terms 
in the last row of Fig. |12[ 

Per definition in ( |3.16| ), M^j^ p,p„f^,f^„ is to be constructed such that its matrix product 
with Mph^piy reproduces Mph^pipnyyi . A low-order manifestation of this is easily verified 
with choosing for M'^plpipi'h'h" uniform limit diagrams shown in the first row of Fig. |T2 



/M-(I) CA 
^^■^ph,p'p"h'h" 



Sh,h'{ph"\T<id\p'p") - Sp,p:{h'h"\r^a\hp") + {p'h') o {p"h" 

+ 5^(p/i'lrddbi/)(pi/ilrdd|V) -$^(p/^lrdd|/iip')(/ii/i'lrdd|V') (B5) 



pi 
1 

N 



fdd(g")(5M'-V) + {p'h')^ip"h") 
+ ^fdd(g")rdd(g') (^h+q' - ?^h+q" 



(B6) 



where the term originating from triplet correlations has not been spelled out. 

Generally, M'^^^i^nyy, is represented by the subset of Mph^pip"h'h" diagrams that can not 
be cut into two pieces, one connected to the labels ph and the other to p'p"h'h'\ by cutting 
either two exchange lines, or cutting the diagram in an internal point. The third row of 



Fig. [12] shows such contributions. 



niihiyi depends non-trivially on three particle and three hole quantum numbers. We 
define the localized version as its Fermi sea average, Eq. ([4.5|) , 



_1 V M^^^""^ 



Sea' 



q,q'+q" 



S{q')S{q" 



Si'\q,q',q") 



[S,iq')S,iq") J S,iq)SM)SAq") S,iq')S,iq 



S{q')S{q") ^ , , , 
+ ^7^^^:hr^us{q,q\q" 



.(B7) 



Here, the relationship (|A4|) was used for the connection between rdd(q') and S{q), and 

Sf^ (g, q', q") = T7 ^h^h-q [^h+q' - ^h+q"] (B8) 



is the three-body static structure function of non-interacting fermions. 
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FIG. 12: Dia grams of Mph,p'p"h'h" in convolution approximation ( p33| ). Graphs obtained by 
exchanging the pairs {p'h') and {p"h") are to be added. The last row shows some diagrams 
containing ground state triplet correlations (shaded triangle), all of these contribute to ^phpipUh'h"- 



3. Three-body vertices 



We now apply the localization procedure ([4.5|) to the three-body vertices. Starting with 



( p.35| ), we have 



^q'q",0 — -"-VV/ 



I y 



hh'h" 



^phpihifl ^^^p'p"h'h",pihi 



pihi 



(B9) 



As discussed in Sec. [IIIB| , the Euler equations (|2.5|) for the ground state optimizations ensure 
that the Fermi sea average (|3.7|) of -f^ppy/i/i'/i'/^o vanishes. For the matrix elements H'^hpiyfl 
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Eqs. O-dD yield 



Gph + <^p'h' — 2 cY^N 



fdd(g)- 



(BIO) 



Therefore, using ( |B6D for M'^^l^^i^n^y, 



jya ^p'p"h'h",0 - jV3 ^phpihifi ^^^p'p"h'h",pihi 

hh'h" hh'h"pihi 



^(?) \ Mil) 



h'h"h 



p'p"h'h",{hi-g)hi 



6. 



,2 c(3)/ 



L5p(g05F(g") 
S^-'iq, q', q") + q • [q" S^q') + q' S,iq")] S^q) 



(Bll) 



This term vanishes when q and g' are larger than 2kp. It is also zero if the matrix element 
Hphpihi,o Eq. (pil| ) is replaced by its Fermi sea average. We therefore expect this term 
to be small, in particular since it has no analog in the Bose limit. Note also that triplet 
ground state correlations do not contribute to this term. Dividing by the normalization 
factors Sp{q) Sp{q') Sp{q") leads to the result ( ^4.13| ). 

To calculate a localized version of the vertex Kph^pip//h'h", Eq- ( |3.34| ), we need 



N S,{q)S,iq')S,{q 



hh'h" 



^ph,p'p"h'h" 2^^ph,pih-t ^^^pihi,p'p"h'h 
Plhi 

(B12) 



with 



Hph,p'h' - '^q,q'|'^/i,?t' eph + ^ ^ph + <^p'h' - 2^gr^ 



fdd(g) . 



(B13) 



We first separate the contribution that survives in the boson limit. Starting with the 
identity 



\ '^p'p"h'h") -F Pq'Pq" - ^ a|,,+^,+g,,a^, (nh'+q" - nh'+q') |$o) 



(B14) 



h'h" 



we have 



-f^p/i,p'p"/i'/i" - (^o I Pq-f^'Pq'Pq"| ^o) - ^ (""-h'+q" " '^h'+q' ) -^^p/i, 



h'+qh' 



(B15) 



hh'h" 



hh' 



Postulating that three-body correlations have been optimized we can simplify the first term 



2N 



[Pq,^'],Pq'Pq" 



2m 



q" S{q') + q' S{q") 



(B16) 
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For the form ( |B13| ), the second term in (pi5| ) is 



- ^ Hph,h>+qh' (r^h'+q" - nh'+q') = — q " q" 5p(g') + q' Sp{q" 



hh' 



2m 



+ ^rdd(g) [g^ Si'^\q, q', q") + q ■ |^q" S'p(g') + q' S'p(g")J S^{q) 

r(i) 



(BIT) 



The remaining term of Kq^q,qn in ( p35| ), -Y.p^h^H'ph,p^h^Mp^hi,p'p"h'h"^ contains contri- 
butions originating from the diagonal and the off-diagonal parts of H'^^^^^^, Eq. ( |B13| ). The 
off-diagonal part is identical to the expression (|B11D , whereas the contribution from the 
diagonal term gives 

S{q')S{q") 



1 /z^ r 

— ep,, M^J = — q- q"5p(g') + q'^F(g" 



hM'h" 



2m 
2m 



[S,{q')SM') 



S{q')S{q")u,{q,q',q"). 



(B18) 



Collecting the individual contributions we obtain Eq. ( ^.121 ). 



4. Four-body coupling matrix element 



In Eq. ( p.24|) we have defined the irreducible four-body coupling matrix element 



^pp'hh' p"p"'h"h"'- Again, "irreducible" means that in the diagrammatic representation left 
and right arguments can not be separated by cutting a particle and a hole line. In analogy 
to the Bose case the "convolution" ( "uniform limit" ) approximation is obtained by retaining 
the leading order diagrams 

^pp'hh',p"p"'h"h"' = ^ph,p"h" Mprh',p"'h"' + Mph,p'"h"' Mp/h',p"h" ■ (B19) 

This contains all diagrams with up to two correlations. A consistent improvement of the 
convolution approximation involves an infinite resummation. For bosons [0 this had only a 
marginal effect. We expect a similarly small improvement for fermions. 

The approximation for Kppihh',p"p"'h"h"' consistent with ( P19|) is to keep all diagrams 
containing only one correlation function rdd(?"), 

^pp'hh',p"p"'h"h"' = ^p,p"^h,h" eph Mp'h',p"'h"' + 5p',p"'Sh',h"' ^p'h' Mph^p"h" 

+ {p"h" ^ p"'h"'} . (B20) 
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Note that both Mp^,^^, p,,^,,,^,,^,,, and K^p^^^^, ^„p,„j^uh,„ contain exphcit particle- and hole-labels. 
Again, we no longer spell out the superscript "CA" in the following. 

A word is in order about the symmetry of both quantities. Eqs. (|B19|) and ( |B20| ) show 
that both operators are the sum of two term that differ from each other merely by the 
interchanging {p"h" p"'h"'}. We have discussed in connection with Eq. (|3.42|) that it is 
legitimate to replace ^pp'/i/i',p"p"'/*"/i"' Kpp'hh',p"p"'h"h'" by their asymmetric form. 



Appendix C: Pair propagator 
1. Pair energy matrix 

A priori, Epp/hh',p"p"'h"h"'{(^) is a function of four hole and four particle momenta as well 
as the energy. In the uniform limit approximation we can, however, express the inverse in 
terms of two-body quantities. From ( [B19[ ) and ( |B2CI| ) we obtain the pair energy matrix 



Eppihh' ,p"p"'h"h"'{^) — {JT^ + iv)^ph,p"h" Mpih',p"'h"' 

— {Sp^p"6h,h" <^ph) Mp/hi^piiih'" — Mph^piiyi {5pi^piii5h',h"' Gp'h') ■ (CI) 
To calculate its inverse, write ( pl|) as 

^^ph,p^hi^^p'h',p2h2^PiV2hih2,p"p"'h"h"'{^) = {^ + 'i'n)^p,p"^h,h" ^p',p"'^h,h"' 

P1/11P2/12 

~ {■^ph,p"h" ^p"h") ^p',p"'^h,h"' — ^p,p"^h,h" (Mpili p,,,^,,, ep'/'h'") (C2) 
Use now, for two commuting operators A, B 



{hw + iri) - A - B 





rdhw' r 


J 


/ 27ri I 



(W+i?7)-A h{uj-uj'+iri) -B 



(C3) 



which can be proved by series expansion. Consequently, we have 

00 

dhw' 



^pp'hh',p"p"'h"h"'i^) 



27ri 



^■ph,p"h" i^') f^p'h',p"'h"' {id — Lu') 



with 



i^ph,p'h'{^) = [{h^ + ir])Mph^pih' - Spp'Shh'eph] ^ ■ 



(C4) 



(C5) 
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For our choice ( ^4.6|) of Mp'h',ph, we can calculate Hph,p'h'i^) analytically, 

l^ph,p'h'[(^) = T — — (^6) 

nw ~ Cph + ir] 

1 ^fdd(g)/iV 1 



where K^{q]u) has been defined in Eq. ( ^.41) . 

According to Eqs. (|3.43|) and (|4.14| ), the dynamic parts of the interactions are obtained 



from matrix products of E^^,^j^i ^,ipn,f^,ii^n,{ijj) as given in (|C4|) with the three-body vertices 
( [4.12| ) and (^4.13|) . The latter being local functions, only sums over the hole states enter 

oo 

E-\q,,q2,u) = ^ J2 E^!p,h,h„p[p',h[h'S'^) = - j ^ <qi;u;') K{q2;u-uj') (C7) 

hih2h[h'2 _^ 

with 

= -T7 2^^ph,p'h'{uj) = -. (C8) 

Using Kramers- Kronig relations, we obtain the useful alternative representation 

CO 

f d{huJi)d{huj2) ^mK{qi;uji)^mK{q2;uJ2) 
J Tc' nwi + nw2 — nw — irj 



2. Properties of the pair propagator 

a. Properties of K{q; oj) 

The structure of n{q] u) resembles that of x{q] ^) iii the RPA. It features a particle-hole 
continuum Kcont(Q'; 1^), and, possibly, a "collective mode" with a dispersion relation given by 

1 + Ko{q; u,{q)) huj,{q) fdd(g) = . (CIO) 

We can therefore write 

3mfi;(g, w) = z{q)iT 6{huj — huc{q)) + '^m K,cont{q]'-^) , 

Z{q) = ^-Jl^ . (Cll) 
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K(g, u) satisfies the following sum rules which we write in the suggestive way 

oo 

S^{q) fd{huj) 





oo 

S^q) fd{huj) 



J^!^^rnK{q;uj) = -S{q) (C12) 



J^!^nLu'^mK{q;u) = -t{q). (C13) 



Eq. ( |C12|) is proved by extending the integration to — oo, noting that Ko{q; u) is real on the 
negative u axis. Since Ko{q;u) has no poles in the upper complex plane, we can evaluate 
the integral along a circle, using the asymptotic expansion 

The proof of Eq. (|C13|) proceeds along the same line, subtracting the asymptotic expansion 
of K{q;uj) beforehand. From Eqs. (|C12|) , ( |C13|) it is clear that the analytic properties of 
S'^{q) K{q; u) / Sp{q) are similar to those of the density-density response function x^^^i^i '^)- 
For bosons, the two functions coincide exactly: Identifying rdd(g') = S{q) — 1 and Sp{q) = l, 
n^{q] oj) consists of a single mode, so that 

"°^^'"^= ;^c.+ir)-t(g) ' ^^''^ 
Figure ^ further confirms this similarity for ^He at saturation density. Expectedly, a 
solution of Eq. ( P10|) is found to lie within a few percent of the RPA zero sound mode. 



h. Properties of E ^{q,q']uj) 

Equations ( P12|) and ( P13|) lead to the sum rules for the pair propagator. 



— SmE = (C16) 

|^r^3,„E-'(,.,» = . (C17) 

— OO 

The proof of (|C16|) is best carried out starting from the representation (|C9D , 



— oo 

oo 



oo oo oo 

/d(hoj) 1 fdhui „ / s f dihw) „ / \ \ 
\smE [qi,q2]U)) = — / ismK{qi;ijJi) / ismK[q2;co — Ui) \L18) 
71 J J TT 
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FIG. 13: Imaginary part of the scaled propagator S (q) K{q,uj)/ Sp{q) (left) and of X {q 
(right) at the density p = 0.0166 A~^. The black squares show, for reference, the Feynman disper- 
sion relation e{q). 

The hw integral in the last term can be extended to — oo since '^mK,{q; u) is real on the 
negative w-axis. 

If Eq. (|C10|) has a solution, the pair propagator has a collective mode. From ( |C11|) we 
obtain 

Q'mE''^{qi,q2;uj) = n z{qi) z{q2) 6{huJc{qi) + ^c{q2) - hu) . (C19) 

This is the origin of two-phonon excitations, or the double-plasmon in charged systems. 

The two-particle-two-hole band consists of three parts which may overlap. The first one 
is the continuum-continuum (c-c) coupling, where the contribution of each /t(g, u) in (|C7|) 
comes from its particle hole band. This defines the two-particle-two-hole "tube" in (g, q'; u) 
space. Its boundaries are 



emin(g) + emin(g') < Hu < emax(g) + emax(g') , 



(C20) 



where Cmin and Cmax denote the upper and lower border of each single-particle-hole band, 
respectively. 

The other two parts of E~^{q^ g'; uj) arise from continuum-mode (c-m) coupling, they are 
identical apart from interchanging q and q' . Their boundaries are 



emin(g) + ^cm(g') < fiW < e^^{q) + hWcmiq') ■ 



(C21) 
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Finally, we consider three limits of the pair propagator. First, in the non-interacting case, 
rdd(Q') = 0, we simply obtain a sum over two-pair energy denominators 

E^\q, q'- uj) = - ! ^f^oiq'; - uj')K^{q; = ^Y] , (C22) 

i.e. the two-particle energy denominator appropriate for perturbation theory in a weakly 
interacting Fermi system. 

Second, ( |C15|) reproduces the energy denominator appearing in the boson theory. 



^ ' " S{q)S{q') nuj + ir^-l{q)-e{q') ' 



Finally, we consider the "collective" or "uniform limit" approximation. Following (|A6|) 
we replace Ko{q;ijj) by that single-pole approximation which ensures its correct and 
sum rules. This gives 



and 

The boson limit as well as the collective approximation demonstrate the effect of cor- 
relations: The single-particle energies get shifted and form a band around the "Feynman- 
spectrum". Note that the collective approximation satisfies the sum rules ( |(J16| )-( |CT7| ) ex- 
actly. 

c. Pair propagator for charged systems 

For charged systems, the dispersion of the solution of Eq. ( |U10|) has, unlike the plasmon, 
a term that is linear in the wave number: 

«-c(<j) = + I f - f f V + 0(9') . (C27) 



6 kp Ahwp ykp 
For the strength of this mode we obtain 



.i.,-.fe)) = ^-3^^. (C28) 
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Hence, to leading order, for the pole of E ^{qi,q2', in ( |C19|) we obtain 



^mE {q , q ; u) = —TT 



l^_lf] s(n^-2nw,-^-^f] ^sq'^0. (C29) 
16 32 A;f / V ^ 3 kp ' 



Note that the location of double-plasmon pole contains, in leading order in the momentum 
transfer, no information on many-body correlations. 

Appendix D: Large momentum limit 



For large momenta, S{q) — 1 falls off at least as q ^. The vertices ([4.12|) and ( [4.13| ) fall 

off as g"^ and as respectively, hence we have 



g'g",o ~ 



5p(g')5p(g") 2m 
0. 



q-q'Xdd(g') + q-q"^dd(g") 



(Dl) 



As a consequence, M4(g;0) is negligible for large momenta, and only the first term in 
Eq ( ^TFD ) contributes to W^{q\^). 

For large q either g' or g" (or both) must be large, (let g">g', the symmetry in q' ^ q" 
just yielding a factor of two). Since Add(g) falls off for large g, the dominant contribution 
of ( pl|) then arises from small g' and we can write 

2 



2 _ 



q-q'Add(g') i5-^(g',g";0) 



Ml 
3 



-\ 2 



7v^dd(g') 



E-^(g',g;0). 



We now use the representation ( |(T^ ) for the pair propagator 



E-i(g',g;0) 



dhjjj' 



TT 



3fJe fi;(g', a;') 3m fi;(g, — cu') 



Since K,^{q^kF] u) = l/{fvjj — t{q) + irj) we have 



Consequently, 



K(g— ^-oo; u) 



E-\q',q^oo;0) 



1 



S'(g) hu! — £:(g) + irj 



^eK{q',-\e{q)) 



Siq) 
t(g) S{q') 



(D2) 



(D3) 



(D4) 



(D5) 
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where the last equahty follows from the high-frequency limit KP{q'; u) — )■ Sp{q')/uj. Insertion 
into dug) yields 

iy^(g^oo,0) = \xUq')\ , (D6) 

q' 



which together with Eq. ( A12 ) gives the result ( 5.14 ). 



Appendix E: Sum rules 

For bosons, the uj^ and sum rules (|1.4|) and ( |1.5|) are satisfied exactly [l^ in the 



sense that the result of the frequency integration is independent of the level at which pair 
fluctuations are treated. This feature provides an unambiguous method to determine the 
static particle-hole interaction V^-J^q) through the sum rule (|1.4|) from the static structure 
function. 

The proof of the nii sum rule is identical to the one for bosons: Due to the symmetry 
we can write 

1 [°° 

mi = 9m / d{hu) hu x{(l'i ■ (El) 

All poles of xiQ'i ^) cire in the lower half plane, allowing to close the integral in the upper 
half plane. For large u we have, however, 

xo{q;uj) -x^^^iq;^) cc^"^ xo(g;c^) - x(g;c^) oc w"^ (E2) 

since 

^ ~ const 

= V^ph(g) H as w oo. (E3) 

We have therefore 

/oo poo roo 

d{huj) hoo xil', = '^fn / (i(^) /za;x^^^(g; cj) = 53m / d{huj) hoo xil', ^) ■ 
■oo J —OD J —OO 

(E4) 



For fermions, the frequency integration in (|1.4|) must be carried out numerically, which 
is best done by Wick rotation along the imaginary axis. The result of the integration is no 
longer rigorously independent of the approximation used for the response function. 



Fig. 14 compares the mo sum rule calculated within the RPA and the pair excitation 



theory. Evidently, the discrepancy is very small. One can understand by comparing with the 
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FIG. 14: Result of the rriQ sum rule for ^He at saturated vapor pressure. The purple dashed line 
shows the FHNC S{q) and the blue short dashed line the result from the pair fluctuation theory; 
the dashed green line shows the difference, magnified by a factor of 30 to make it visible. 

boson theory: If we restricted the fluctuation operators 6u^p^ (t) and Su^pp/hh'i^) be functions 
of momentum transfers q = p — h and q' = p' — h', we would end up with a density- density 
response function that is formally identical to that of bosons and would, hence, lead to an 
S{q) that is independent of the treatment of the pair fluctuations. The expectation that 
the inclusion of the particle-hole structure of the two-pair energy denominator makes only a 
small difference is verified in Fig. [1^. Thus, it is also legitimate in the pair-excitation theory 
to obtain the static particle-hole interaction Vp_h(g) from the static structure function S{q) 
through Eqs. and ( [Ol) . 



Appendix F: Implementation Recipe 

This section provides, for the convenience of the reader and easy further reference, a 
compilation of all necessary ingredients to implement the theory. Mostly a summary of 
sections |I^ and |V A| , we deliberately refrain from any explanation to avoid redundancy and 



keep it as compact as possible. 

We have shown in our applications to '^He and the electron liquid that for practical 
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purposes, only one of the local three-body vertices is necessary: 

S{q')Siq") 



^q,q'q" — ' 



2m S,{q)S,{q')S,{q") 



q-q'Xdd(g') + q-q"^dd(g")-r^3(g,g',g") , (Fi) 



where u^^q, q' , q") is the three-body ground state correlation p8[. The effective interaction 
W^{q, u) is then 



W,{q;u) 



2N 5^ '^"'^'^ 



?E-\q',q"-u) 



(F2) 



whereas W^{q,u}) vanishes. Consequently, the components of the (energy-dependent) inter- 
action matrix Vp_h(i^) are 



(F3) 
(F4) 



witha±^[^p(g)±^(g)]/2^(g). 

Finally we need the pair propagator: 

E-^{qi,q2]Uj) = - 

K,{q;u) 



K{qi;uj') n{q2\UJ-u') 

2m 

Ko{q;uj) 



1 + huTM{q)i^o{q;(^) 



(F5) 
(F6) 



with the partial Lindhard functions: 



1 n-p^h 
KJq;uj = — > -T — 



The simplifications of the interactions do not significantly simplify the form 
density- density response function. 



(F7) 
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